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OVERVIEW 



24 K 


D SUCH DATA COULD BE USED TO CREATE AN “EXPERT METEOROLOGIST USING 
EITHER Al OR OTHER TECHNIQUES THAT COULD BE USED IN EITHER REACTIVE OR 
LOOK-AHEAD SYSTEMS TO VARY SENSITIVITY THRESHOLDS AND COORDINATE THE 
INPUTS PROM DIFFERENT DETECTING SYSTEMS. 




OVERVIEW 
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The 

Numerical Field Experiment 

Methodology 


Use a Numerical Model to Simulate a Large 
Population of Physically Plausiable 
Scenarios Similar to What Might be Anticipated 
During a Field Observing Program 
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Goals of Study 


• Better understand how various physical 
parameters interact to determine microburst 

intensity 


Determine probabilities of microburst occurrence 
under a variety of conditions 


• Validate results against observations 



Experiment Design 


• Axisymmetric Numerical Cloud Model 

• Warm Rain Microphysics 

• No Ambient Wind 

• Zero Ambient Humidity 

• Simulate Only the Sub-Cloud Region 

• Continuous Influx of Rainwater at Model Top 
(Cloud Base) 
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Vertical Coordinate 




Radial Coordinate 
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Parameter Space 


Cloud Base Height (0.5, 1, 2, 3, 4 km AGL) 
Sub-Cloud Lapse Rate (70, 80, 90, 100% D.A.) 
Surface Temperature (55, 65, 75, 85, 95, 105 F) 
Cloud Base Reflectivity (20, 30, 40, 50, 60 dBz) 
Rainshaft Radius (0.5, 1.0, 2.0 km) 

All Combinations -» 1800 Simulations 

• Each Run for 20 Minutes 

• Max's & Min's of all Fields Saved Every 10 s 



Model Validation 


• NCSA Bakeoff - Tested Against Some 15 Codes 

• Independent Tests with Krueger Axisymmetric 

and Klemp & Wilhelmson 3-D 



Umax 

Wmin 

Droegemeier 

46.7 

-35.0 

Krueger 

45.1 

-34.8 

Droegemeier 

52.5 

-37.9 

K & W 3-D 

47.4 

-37.9 



(m r *) 








5.0 lO * 2 


4.0 10* 2 


^ 3.0 10* 2 

§ 

5 2.0 10- 2 


1.0 10* 2 
0.0 10 ° 
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Number of Occurrences 



Number of Occurrences Number of Occarreaces 
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Energy (m 2 s* 2 ) 
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Radial Wind Speed (m s - ') 



Temperature Deficit (K) 
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Reflectivity (dBz) 











Lap«e Rate (% of Dry Adloballc) Cloud Bate Hctfbl (k») 







Table S. Slope of linear fit between radial and downdraft wind speeds for 
ail 1800 simulations as a function of rainshaft radius (R) and cloud base height (H), the ratio 
of which is defined as the aspect ratio (R/H). 





Rainshaft Radius R (km) 

Cloud Base H (km) 

Aspect Ratio (R/H) Slope 


2.0 

0.5 

4.00 

-0.60 

1.0 

0.5 

2.00 

-0.78 

2.0 

1.0 

2.00 

-0.75 

2.0 

2.0 

1.00 

-0.90 

1.0 

1.0 

1.00 

-0.93 

0.5 

0.5 

1.00 

-0.97 

2.0 

3.0 

0.67 

-0.95 

2.0 

4.0 

0.50 

-0.97 

1.0 

2.0 

0.50 

-1.03 

0.5 

1.0 

0.50 

-1.10 

1.0 

3.0 

0.33 

-1.04 

1.0 

4.0 

0.25 

-1.06 

0.5 

2.0 

0.25 

-1.16 

0.5 

3.0 

0.17 

-1.17 

0.5 

4.0 

0.13 

-1.17 



Table 6. Probabilities (%) of radial wind speed classified according to model input parameters for 

the set of 1800 simulations. 

aft, indicates that the associated parameter varies among all values used in the Subset 


Input Variable 


LR Tsfc Refl CB Rad 


Wind Speed 


> 10.3 ms' 1 >12.9 ms* 1 >15.4 ms* 1 > 18.0 ms* 1 
(> 20 kts) (£ 25 kts) (> 30 kts) (2> 35 kts) 


ALL ALL ALL ALL ALL 


44.50 35.17 28.22 22.61 


100 ALL ALL ALL ALL 59.11 

90 ALL ALL ALL ALL 45.33 

80 ALL ALL ALL ALL 39.56 

70 ALL ALL ALL ALL 34.00 


45.56 

37.11 

31.11 
26.89 


36.89 29.11 

29.78 24.00 

24.44 19.78 

21.78 17.56 


ALL 

105 

ALL 

ALL 

95 

ALL 

ALL 

85 

ALL 

ALL 

75 

ALL 

ALL 

65 

ALL 

ALL 

55 

ALL 


ALL 

ALL 

48. 

ALL 

ALL 

48. 

ALL 

ALL 

46. 

ALL 

ALL 

45. 

ALL 

ALL 

42. 

ALL 

ALL 

36. 


67 40.67 

00 39.33 

33 37.00 

33 34.33 

33 31.67 

33 28.00 


33.33 26.33 

31.33 26.00 

29.67 25.33 

28.00 23.00 

25.33 20.00 

21.67 15.00 


ALL 

ALL 

60 

ALL 

ALL 

50 

ALL 

ALL 

40 

ALL 

ALL 

30 

ALL 

ALL 

20 


ALL 

ALL 

98. 

ALL 

ALL 

74. 

ALL 

ALL 

36. 

ALL 

ALL 

11. 

ALL 

ALL 

1 . 


89 85.28 

17 62.50 

39 22.50 

1 1 5.28 

94 0.28 


78.61 70.56 

46.94 35.00 

13.61 7.22 

1.94 0.28 

0.00 0.00 


ALL ALL ALL 4.0 ALL 

ALL ALL ALL 3.0 ALL 

ALL ALL ALL 2.0 ALL 

ALL ALL ALL 1.0 ALL 

All ALL ALL 0.5 ALL 


49.44 

53.06 

56.11 

45.00 

18.89 


45.56 

46.39 

46.94 

31.67 

5.28 


39.72 

41.11 

38.33 

21.94 

0.00 


34.44 

35.56 

31.67 

11.39 

0.00 


AIL ALL ALL ALL 2.0 

ALL ALL ALL ALL 1.0 

ALL ALL ALL ALL 0.5 


52.00 43.33 

45.17 37.33 

36.33 24.83 


35.50 29.00 

29.33 25.17 

19.83 13.67 


LR = Lapse Rate (% of Dry Adiabatic) 
Tsfc * Surface Temperature (F) 

Refl = Reflectivity Factor (dBz) 

CB = Cloud Base Height (km) 

Rad = Rainshaft Radius (km) 
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m r 1 ) 




m m § 


Table 7 Probabilities (%) of radial wind speed classified according to model input parameters 

a subset of 1 14 simulations defined below. 

ALL indicates that the associated parameter varies among all values used in the Subset 


for 


Input Variable Wind s P ccd 


Tsfc 

Refl 

CB 

Rad 

> 10.3 m s'l 
(> 20 kts) 

> 12.9 ms'* 
(2 25 kts) 

> 15.4 m s*l 
(> 30 kts) 

> 18.0 m s 
(2 35 kts) 

ALL 

ALL 

ALL 

ALL 

60.53 

41.23 

32.46 

21.05 

ALL 

ALL 

ALL 

ALL 

60.53 

41.23 

32.46 

21.05 

105 

95 

85 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

60.53 

60.53 

60.53 

42.11 

42.11 

39.47 

34.21 

31.58 

31.58 

21.05 

21.05 

21.05 

ALL 

ALL 

ALL 

ALL 

ALL 

60 

50 

40 

30 

20 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

100.00 

66.67 

55.56 

28.57 
25.00 

85.19 

55.56 

22.22 

14.29 

0.00 

66.67 

48.15 

22.22 

0.00 

0.00 

55.56 

22.22 

11.11 

0.00 

0.00 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

2.0 

1.0 

0.5 

ALL 

ALL 

ALL 

85.71 

61.54 

27.27 

64.29 

38.46 

15.15 

57.14 

33.33 

0.00 

42.86 

15.38 

0.00 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

2.0 

1.0 

0.5 

60.00 

66.67 

54.55 

46.67 

47.22 

27.27 

33.33 

36.11 

27.27 

26.67 

25.00 

9.09 


Lapse Rate * 100% Dry Adiabatic 
85 F £ Surface Temperature £ 105 F 
20 dBz £ Reflectivity £ 60 dBz 
0.5 km £ Goud Base Height £ 2.0 km 
0.5 km £ Rainshaft Radius <. 2.0 km 

Urad ^ 5 m s ^ 
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Table I. C o effic ients of multivariate linear regression equations for the radial wind speed (ms* 1 ) 
for ail 1800 simulations ( 1799 degrees of freedom), using various combinations of predictors. 


Shown are the predictors, the estimated standard deviation of model error (m s' 1 ), the leading 
constant and the five regression coefficients. The units of the regressors are as shown in Table 3, 
except for the lapse rate, which is given by its decimal equivalent 


Predictor(s) 

% 

Variance 

Explained 

Est Std. 
Dev. of 
Model Error 


Regression Coefficients 



Const 

LR 

Tsfc 

dBz 

CBHt 

Rad 

LR 

3.7 

9.8 

-3.40 

17.31 

— 

— 

— 

— 

Tsfc 

1.7 

9.9 

5.15 

— 

0.08 

— 

— 

— 

dBz 

60.4 

6.3 

-10.69 

— 

— 

0.55 

— 

— 

CBHI 

3.4 

9.7 

7.50 

— 

— 

— 

1.82 

— 

Rad 

4.6 

9.8 

7.30 


— 

— 

— 

3.43 

dBz, LR 

64.1 

6.0 

-25.38 

17.31 

— 

0.55 

— 


dBz, Tsfc 

62.1 

6.2 

-16.83 

— 

0.08 

0.55 

— 

— 

dBs,CBHt 

65.8 

5.9 

-14.51 

— 

— 

055 

1.82 

— 

dBz, Rad 

63.0 

5.9 

-14.69 

— 

— 

0.55 

— 

3.43 

LR, Tsfc 

5.5 

9.7 

-9.56 

17.31 1 

0.08 

— 

1.82 

— 

LR, CBHt 

9.2 

9.5 

-7.22 

17.31 

— 

— 

— 

LR, Rad 

8.3 

9.6 

-7.41 

17.31 

— 

— 

— 

3.43 

CBHt, Tsfc 

7.1 

9.7 

1.33 

— 

0.08 

— 

1.82 

— 

CBHt, Rad 

10.0 

9.5 

3.48 

— 

— 

— 

1.82 

3.43 

Tsfc, Rad 

6.3 

9.7 

1.14 

— 

0.08 

— 

— 

3.43 

LR. Tsfc, CBHt 

10.9 

9.5 

-13.38 

17.31 

0.08 

— 

1.82 

— 

LR. Tsfc, dBz 

65.9 

5.9 

-31.57 

17.31 

0.08 

0.55 

— 

— 

LR, Tsfc, Rad 

10.1 



-13.57 

17.31 

0.08 

— 


3.43 

LR, CBHt. dBz 

69.6 

5.5 

-29.22 

17.31 

— 

0.55 

1.82 

— 

LR, CBHt, Rad 

13.7 

9.3 

-11.23 

17.31 

— 

— 

1.82 

3.43 

LR, dBz, Rad 

68.7 

5.6 

-29.41 

17.31 

— 

0.55 

— 

3.43 

LR, Tsfc, CBHt, Rad 

15.5 

9.2 

-17.39 

17.31 

0.08 

— 

1.82 

3.43 

LR, Tsfc, CBHt, dBz 

71.3 

5.4 

-35.39 

17.31 

0.08 

0.55 

1.82 

3.43 

LR, Tsfc, Rad, dBz 

70.5 

5.5 

-35.57 

17.31 

0.08 

0.55 

— 

LR, CBHt, Rad, dBz 

74.1 

5.1 

-33.23 

17.31 


0.55 

1.82 

3.43 

LR, Tsfc, CBHt, 
dBz. Rad 

75.9 

4.9 

-39.40 

17.31 

0.08 

0.55 

1.82 

3.43 
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Table 9. As in Table 8, but assuming that the cloud base height is known (359 degrees of 
freedom) and takes the values shown in the left column. 


I % Est Std. 

Cloud Base I Variance Dev. of 

Explained Model Error 


Regression Coefficients 


Const 


LR Tsfc dBz CBHt Rad 


.23 — 1.27 



























Tabic 10. As in Table 8, but assuming that the reflectivity factor is known (359 degrees of 
freedom) and takes the values shown in the left column. 


Radar 

Reflectivity 

% 

Variance 

Explained 

EsL Std. 
Dev. of 
Model Error 

Regression Coefficients 

Const 

LR 

Tsfc 

dBz 

CBHt 

Rad 


60.2 

1.5 

-9.08 


\ ] 1 

— 

•0.46 

1.94 

20 dBz 

24.8 

2.0 

-6.39 

. 1 

1 

— 

— 

— 


1.7 

2.3 

0.98 

I 

0.02 

— 

— 

— 


6.6 

2.3 

3.37 

I 

— 

— 

-0.46 

— 


27.2 

2.0 

0.13 


— 

— 

— 

1.94 


56.3 

2.5 

-18.04 

20.73 

0.03 



•0.01 

2.49 

30 dBz 

37.6 

3.0 

-12.76 


— 

— 

— 



1.8 

3.8 

2.30 

— 

0.03 

— 

— 



0.0 

3.8 

4.89 

— 

— 

— 

-0.01 

IPia 


16.8 

3.4 

1.96 

- 

— 

— 

— 

— 

2.49 


56.5 



3.7 

L 

-25.22 

26.84 

El 



0.67 

3.85 

40 dBz 

30.1 

4.6 

-13.77 

26.84 

lifi 

— 

— 

— 


4.7 

5.4 

3.49 

— 

0.07 

— 

— 

— 


2.5 

5.4 

7.62 

— 

— 

— 

0.67 

— 


19.2 

4.9 

4.56 

; 

— 



— 

— 

3.85 


66.7 

4.0 

•21.52 





3.04 

4.36 

50 dBz 

9.6 

6.5 

-0.27 


| . 

— 

— 



9.3 

6.5 

6.13 


0.12 

— 

— 



32.2 

5.6 

9.51 

I I 

— 

— 

3.04 



15.7 

6.3 

10.81 

WSM 

— 

— 


4.36 


'■mu-.- 

83.2 

3.8 

B39 

9.65 

0.14 



5.85 

4.55 

60 dBz 

1.4 

9.3 

IB 1 

9.65 

— 

— 

— 

— 


7.2 

9.0 

1BE1 

— 

0.14 

— 

— 

— 


65.3 

5.5 

IBB 

— 

— 

— 

5.85 

— 


9.3 

8.9 

19.07 

— 

— 

— 

— 

4.55 
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Tabic 12. As in Table 8, but assuming that the rainshaft radius is known (599 degrees of freedom) 

and taki»s the values shown in the left column. 


Rainshaft 

Radius 

% 

Variance 

Explained 

EsLStd. 
Dev. of 
Model Error 

Regression Coefficients 

Const 

LR 

Tsfc 

dBz 

CBHt 

Rad 


77.9 

3.8 

•29.05 

15.27 

0.05 

0.47 

1.02 

_ 

0.5 km 

4.6 

7.8 

- 4.60 

15.27 


— 

— 

— 


0.9 

7.9 

4.86 

— 

0.05 

— 

— 

— 


69.7 

4.4 

- 10.41 

— 

— 

0.47 

— 

— 


2.7 

7.9 

6.25 

— 

— 

— 

1.02 

— 


77.1 

4.9 

-37.28 

18J3 

0.07 

0.59 

1.91 



1.0 km 

4.0 

10.0 

- 3.88 

18.33 

— 

— 

— 

— 


1.5 

10.2 

5.81 

— 

0.07 

— 

— 

— 


65.9 

6.0 

- 11.80 

— 

— 

0.59 

— 

— 


5.7 

9.9 

7.69 

^ 

— 

i 

1.91 

— 


75.7 

5.4 

•39.82 

18.35 

0.11 

0.59 

2.53 

— 

2.0 km 

3.6 

10.7 

- 1.73 

18.35 

— 

— 

— 

— 


3.2 

10.7 

4.78 

— 

0.11 

— 

— 

— 


59.9 

6.8 

- 9.85 

— 

— 

0.59 

— 

— 


8.9 

10.3 

8.55 

— 

— 

— 

2.53 

— 






















































Table 13. As in Table 8, but assuming that the lapse rate is known (449 degrees of freedom) and 

takes the values shown in the left column. 
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Preliminary Comparison 
with JAWS Observations 


Model Data Subset Satisfying the Following: 

U> lOms-1 

2.0 km < Cloud Base < 4.0 km 
75 F < Sfc Temp ^ 95 F 
20 dBz < Reflectivity < 40 dBz 
90% D. A. < Lapse Rate < 100% D.A. 
0.5 km < Rainshaft Radius < 2.0 km 


72 Simulations (vs 186 JAWS events) 



Table 6. Probabilities of horizontal wind speed computed from the model data and JAWS 

observations. 


Wind Speed 

Dataset 

> 20 kts 

> 25 kts 

> 30 kts 

>35 kts 

JAWS Observations 

98.39 

76.88 

42.47 

23.12 

All 1800 Simulations 

44.50 

35.17 

28.22 

22.61 

tModel JAWS Subset 

98.61 

66.67 

40.28 

25.00 


t Criteria: radial wind ^ 10 m s ^ 

2.0 £ cloud base height £ 4.0 km 

75 F £ surface temperature £ 95 F 

20 dBz £ cloud base reflectivity £ 40 dBz 

90% dry adiabatic $ lapse rate £ 100% dry adiabatic 

0.5 km £ rainshaft radius £ 2 km 

(72 cases) 
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Preliminary Conclusions 


• The Model Solution Space Contains Considerable 
Variability 

• Solution Behavior is Physically Consistent 

• Reflectivity is the Dominant Influence Most of 
the Time in Determining Radial Wind Speed 

• Gross Statistical Comparisons with JAWS 
are Encouraging 

• The "Numerical Field Experiment" Approach 
Appears Well-Suited to this Problem 

• Several Avenues for Practical Application 
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Ongoing Work 


• Inclusion of Ice Processes and Shallow Stable 
Layers 


• "Pointwise" and Overall Statistical Comparisons 
with JAWS and Orlando TDWR OT&E Data 


• Refinement of Parameter Space 


• Evaluation of Probability Calculations for 
Operational Forecasting 

• Further Development of Parametric Relationships 
Between Model Outputs and Inputs 
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A "Numerical Filed Experiment" Approach for 
Determining Probabilities of Microburst Intensity 

Questions and Answers 

Q: Kim Elmore (NCAR) - 1 am fascinated by the work you did, and specifically how much of a 
predictor reflectivity was in general. For the equivalent JAWS simulations, how good of a 
predictor was reflectivity for intensity?. Do you remember? 

^"?oT 8emeier (Universi,y of 0k,ahoma ) - We haven’t actually done that break down 
for the JAWS data yet. Based on the other results, I would say it was probably a very strong 
influence. . 

Kim Elmore (NCAR) - Our experience in JAWS was that it wasn’t a veiy good predictor of the 
outflow we would see. That was one of the major conclusions. Reflectivity, for the JAWS data 
was not a good predictor. * 

Kdvin Droegemeier (University of Oklahoma) - 1 just had a student finish a Master’s Thesis on 
a study of Orlando cases. What we have found and what tended to make that conclusion seem 
plausible in light of the fact that the reflectivity factor was the same and apparently similar 
environments, was the fact that low level effects like low level inversions in stable air change the 
outflow intensity. It does not take much stabilization in low levels to really change the outflow 
intensity. In fact you will see that if you stay around and look at the animation sequence. Now I 
know Fred and others of us who run models have actually dropped globs of rain into stable air for 
a long time. What this student did, for the first time, with a 3-D cloud simulation with ice, 
showed that the storms themselves forming beneath the low level stable air, were virtually’ 
unaffected by it. Once the rain came down and the outflow hit, that is when the radial winds were 
really diminished by virtue of the stable air. 

Kim Elmore (NCAR) - We have seen, when gust fronts go by, that often once you stabilize a 
relatively deep part of the boundary layer that they do not make it through any more. 

Kelvin Droegemeier (University of Oklahoma) - Yes; it doesn't have to be but maybe four or 
five hundred meters. It can be pretty shallow. One of the limitations of our study is that we did 
not put in the shallow stable layer. If you consider the parameters, you have to figure how thick 
is a layer and how cold is it. That is two more parameters and that would run it up to 20 000 
simulations. 

Kim Elmore (NCAR) - In JAWS, we also did not stratify the cases as to was there cold surface 
layer or not. 

Kelvin Droegemeier (University of Oklahoma) - That is a tough thing to do. Usually they just 
skim the mesonet and you just never know. 

Kim Elmore (NCAR) - Our gut feeling was; once we had a cold surface layer, if it had been 
around long enough to deepen, however much it had to deepen, and we did not even really know 
how much that was, that tended to shut it off. 
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Kelvin Droegemeier (University of Oklahoma) - 1 completely agree, I think that is the 
controlling influence for attenuating that stuff. 

Q: Brae Bracalente (NASA Langley) - The reflectivity values you referred to are at the cloud 
level right at the high altitude where the rain first starts to fall. Did you do any correlating with 
the reflectivity levels at the outflow region? We have found that the peak shear did not 
necessarily occur where the heaviest rain was, 

A: Kelvin Droegemeier (University of Oklahoma) - No, we have not done that but we could. 
In fact, if you stay around for the video tape you will see that the reflectivity near the ground or 
not to far above the ground is less than it is at cloud base height. I should have mentioned, in the 
absence of having ice in this case, we are assuming that once the precipitation falls below cloud 
base the only stuff that is important for the forcing that occurs that drives that microburst is the 
stuff that happens below cloud base. Obviously, that is not always the case, but that was the 
assumption in this case here. The reason we did that is so we do not have to consider all the 
possible soundings and wind profiles and everything that happens above cloud base where there is 
a lot of variability. So that was the other assumption. 

Q: Fred Proctor (NASA Langley) - Define what you mean by your cloud base height in your 
model studies? 

A: Kelvin Droegemeier (University of Oklahoma) - It is basically just the height at which the 
rain begins to fall, that is the simplest explanation. 

Q: Fred Proctor (NASA Langley) - Is that the top of your model then? 

A: Kelvin Droegemeier (University of Oklahoma) - Yes, and the vertical velocity is zero 
there. 
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